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Viscous hydrodynamics is commonly used to model the evolution of the matter created in an 
ultra-relativistic heavy-ion collision. It provides a good description of transverse momentum spectra 
and anisotropic flow. These observables, however, cannot be consistently derived using viscous 
hydrodynamics alone, because they depend on the microscopic interactions at freeze-out. We derive 
the ideal hydrodynamic limit and the first-order viscous correction to anisotropic flow (u 2 , wa and 114) 
and momentum spectrum using a transport calculation. The linear response coefficient to the initial 
anisotropy, Vn{pT)/e-n., depends little on n in the ideal hydrodynamic limit. The viscous correction 
to the spectrum depends not only on the differential cross section, but also on the initial momentum 
distribution. This dependence is not captured by standard second-order viscous hydrodynamics. 

The viscous correction to anisotropic flow increases with px, but this increase is slower than usually 
assumed in viscous hydrodynamic calculations. In particular, it is too slow to explain the observed 
maximum of Vn at pr ~ 3 GeV/c. 


I. INTRODUCTION 

Relativistic viscous hydrodynamics [10 is the state 
of the art for describing the evolution of the strongly- 
coupled system formed in an ultrarelativistic heavy-ion 
collision at RHIG or LHC. It has long been realized 0 
that ideal hydrodynamics naturally explains the large 
magnitude of elliptic flow @, 0 . However, the system 
formed in such a collision is so small that deviations 
from local thermal equilibrium are sizable, resulting in 
the inclusion of viscosity 0 in hydrodynamic calcula¬ 
tions. Viscosity twically reduces the magnitude of ellip¬ 
tic flow by 20% 0. Viscous effects on higher harmonics 
of anisotropic flow 0 0 , such as triangular flow 0 , are 
even larger (TTj |. 

Even though there is a consensus that viscosity mat¬ 
ters, the calculation of viscous corrections to observables 
is not yet under control. The reason is that viscosity af¬ 
fects not only the space-time history of the fluid [ 10 , but 
also the momentum distribution of particles at “freeze- 
out” , which has an off-equilibrium part proportional to 
viscosity [H, [T3|. Viscous hydrodynamics itself does not 
fully specify this off-equilibrium part. The only require¬ 
ment is that the system of particles should generate the 
same ener gy- momentum tensor as the fluid just before 
freeze out [l5|. This requirement, however, does not con¬ 
strain the dependence of the relative deviation to equilib¬ 
rium on the momentum p in the rest frame of the fluid, 
which is essentially a free function xip)- This function is 
not universal, and involves the differential cross sections 
between constituents [l0. It is typically put by hand in 
hydrodynamic calculations. 

The common lore is that effects of viscosity are 
more iimortant for particles with larger transverse mo¬ 
menta [0. This is due to the fact that mos t hy drody- 
namic calculations use the “quadratic” ansatz [Q x{p) ^ 


p^. While this choice generally results in a improved de¬ 
scription of experimental data 0 (see howeve r fisl ). it 
is not supported by any theoretical argument |16l |. We 
evaluate viscous corrections to observables (specifically, 
transverse momentum spectra and anisotropic flow) by 
solving numerically a relativistic Boltzmann equation. 
We simulate relativistic particles undergoing 2 —>• 2 elas¬ 
tic collisions with a total cross section crtot- In the limit 
Ctot —>■ + 00 , a generic observable /(cr) can be expanded 
in powers of 1 /crtot: 

/Kot) « + —+ . (1) 

The leading term is the limit of infinite cross sec¬ 
tion, which corresponds to ideal hydrodynamics in the 
limit of a vanishing freeze-out temperature [13, [ 3 . The 
next-to-leading term 6f is a viscous correction: since the 
shear viscosity 77 scales like I/ctot 3 , this correction is 
proportional to ryO 

We evaluate and Sf by solving numerically the rel¬ 
ativistic Boltzmann equation for several large values of 
the cross section crtot- Our primary goal is to illustrate 
by an explicit calculation how the viscous correction to 
anisotropic flow depends on transverse momentum px, 
and to what extent this dependence is sensitive to the 
structure of the differential cross section. We do not 
mean here to carry out a full realistic simulation of a 
heavy-ion collision. In particular, for sake of simplicity, 
our transport calculation uses massless particles which 
supply the possibility of having only shear viscosity with 


^ Note that in hydrodynamic calculations, 5f often denotes the vis¬ 
cous correction at freeze-out. Here, Sf is the full viscous correc¬ 
tion, which also contains a contribution from the hydrodynamic 
evolution. 
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no bulk viscosity (20l - [^ . The resulting equation of state 
is harder than that of QCD near the deconfinement 
crossover |24| . This results in larger Vn and harder pT 
spectra. 

We also study the dependence of observables on ini¬ 
tial conditions. In second-order viscous hydrodynamics, 
the evolution is completely specified by the initial value 
of the energy-momentum tensor |I 3 | • At the mi¬ 

croscopic level, however, the initial momentum distribu¬ 
tion contains more information than just . Usual 
viscous hydrodynamics assumes that this additional in¬ 
formation is washed out by the system evolution. Our 
simulation provides a means of testing this assumption, 
by constructing two different initial conditions with the 
exact same , and comparing the observables at the 
end of the evolution. 


II. INITIAL CONDITIONS AND EVOLUTION 

Initial conditions follow Bjorken’s boost-invariant pre¬ 
scription [^ . but with a finite extent in space-time ra¬ 
pidity —2.5 < rj < 2.5. The evolution is started at time 
To = 0.6 fm/c [i^ after the collision. The initial con¬ 
ditions of the Boltzmann equation are specified by the 
one-body density f(x,p) in coordinate (x) and momen¬ 
tum (jp) space at time tq. 

The initial density profile in transverse coordinate 
space (x, y) is taken from an optical Glauber calcula¬ 
tion for a central Au-Au collision at i/s = 200 GeV, cor¬ 
responding to the top RHIC energy. This initial density 
is azimuthally symmetric, so that anisotropic flow van¬ 
ishes by construction. We introduce anisotropy artifically 
by deforming the initial distribution, thus mimicking an 
initial state fluctuation M- In hydrodynamicSjOne typi¬ 
cally deforms the initial energy density profile [III, . In 
a transport calculation, where the initial conditions are 
specified by the initial positions of particles, it is simpler 
to just shift these positions by a small amount. Introduc¬ 
ing the complex notation z = x + iy, in order to generate 
flow in harmonic n, we shift z according to 

z^z + 6z = z — ( 2 ) 

where z = x — iy, and a is a real positive quantity chosen 
in such a way that the correction is small. This transfor¬ 
mation is invariant under the change ^ —>■ i.e., 

it has 27r/n symmetry. Therefore, to leading order in ix 
the only nonvanishing anisotropic flow coefficient is Vnu 
The initial eccentricity in harmonic n is defined for n > 2 


^ Note that the deformation used in Refs. IllI |2^ is singular at 
the origin z = 0 for n 2. By contrast, the deformation defined 
by Eq. 1(2} is regular at the origin, but has a singularity at large 
\z\ (see also [3). On the x axis, the singularity is located at the 
point where dSxfdx = —1. If a. is small, this singularity occurs 
at a point where the density is low. 


by [21,133 


J2jizj + Szj) 

\Zj+SZj\ 


(3) 


where the sum runs over all particles with initial position 
Zj. Inserting Eq. m into Eq. and using the fact 
that the distribution of Zj is azimuthally symmetric, one 
obtains, to leading order in a and for a large number of 
particles 


E, \Zj\^ 

/^2(n-l)\ 

~ na -—;—;—, 


= na- 




Ej lul 


(4) 


where, in the right-hand side, angular brackets denote 
an average taken with the optical Glauber profile, and 
r = l^l = Y^x^ + 2 /^. We carry out simulations for 
n = 2,3,4. For each n, we fix a in such a way that 
En = 0.2. We have checked that this value is sufficiently 
small that the response is linear m- In hydrodynamics 
with fluctuating initial conditions, V 2 and are deter¬ 
mined to a good approximation by linear response to 
the eccentricity in the corresponding harmonic (El . . 
On the other hand, X 4 is the superposition of a linear 
term (33j |. and a nonlinear term induced by V 2 [^ . The 
present study only addresses the linear term. 

In momentum space, we consider two different types of 
initial conditions. The first case is a thermal Boltzmann 
distribution 


dN/(fp(xexp{—p/T). (5) 

In addition, one requires that the temperature T and the 
local density n be such that n/T^ is a constant through¬ 
out the transverse plane, as in a thermal gas of massless 
particles with zero chemical potential. Thus this initial 
condition is that of a hydrodynamic calculation with zero 
chemical potential (23|. The initial temperature at the 
center of the fireball is Tq = 340 MeV. The second case 
is a constant distribution: 

dN/(fp oc 6{po — p)- (6) 

The maximum momentum pq and the proportionality 
constant are chosen such that the particle density n and 
the energy density e are the same as with the previous 
initial condition. Similar initial conditions have been pre¬ 
viously used in transport calculations in order to mimic 
the effect of saturation in high-energy QGD [3l,|33. Both 
types of initial conditions have exactly the same energy- 
momentum tensor. 

The evolution of the system is determined by the rela¬ 
tivistic classical Boltzmann equation. We use a relativis¬ 
tic transport code developed to study heavy-ion collisions 
at RHIG and LHG energies [13 [tH 441|| , which uses the 
test-particle method. The collision integral is solved by 
using Monte Garlo methods based on the stochastic in¬ 
terpretation of transition amplitudes [l 3 , [13 Ell- The 
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total cross section is fixed throughout the evolution. The 
differential cross section is 


da 1 

oc - 

dt (^t — mjj) 


(7) 


where t is the usual Mandelstam variable. This differen¬ 
tial cross section is typically used in parton cascade ap¬ 
proaches [13, m, l43 - l44| and by symmetry the w-channel 
is included. The limit tou —?► oo corresponds to an 
isotropic cross section. The opposite limit, where mo 
is smaller than the typical particle energy, corresponds 
to a forward-peaked cross section. 

In a transport approach, one can follow the eyolution of 
the system until the last collision, but this is numerically 
expensiye. Instead, we choose to follow the system until 
a fixed final time, and check stability of our results with 
respect to this final time. The calculations presented in 
this paper are carried out with a final time tf = 8 fm/c, 
but we haye checked that the momentum spectra are un¬ 
changed if we extend the final time to tf = 12 fm/c. Since 
the initial density profile possesses 4> ^ ~ 4’ symmetry 
(where is the azimuthal angle), we define anisotropic 
flow as Vn = (cos n(j )), where angular brackets denote an 
ayerage oyer particles at the end of the eyolution. 

Throughout this paper, we carry out four sets of calcu¬ 
lations: Three sets with the thermal initial distribution 
dS]) and the yalues mo = 0.3 GeV, mo = 0.7 GeV and 
an isotropic cross section, and a fourth set with the con¬ 
stant initial distribution m and an isotropic cross sec¬ 
tion. Thus we study how results depend on the initial 
distribution and on the differential cross section. 

For each set of parameters, we perform different cal¬ 
culations for the following set of total cross sections: 
CTtot = 20, 25, 30 and 35 mb. Results for anisotropic 
flow Vn [H are shown in Fig. [1] The dependence of these 
obseryables on l/utot is essentially linear, corresponding 
to the regime where yiscous hydrodynamics applies. In 
order to improye the accuracy, we fit these results with 
a polynomial of order 2 in l/utot and extract the ideal 
hydrodynamic limit and the first yiscous correction us¬ 
ing Eq. (HD: specifically, the ideal hydrodynamic limit 
is the extrapolation to 1/crtot —t 0 and the yiscous 
correction Sf is the slope at the origin. 


III. IDEAL HYDRODYNAMICS 

We first study the ideal hydrodynamic limit, defined by 
the extrapolation atot -l-oo in the Boltzmann equation. 
Ideal hydrodynamics corresponds to local thermal equi¬ 
librium [^ . One expects thermalization to wash out de¬ 
tails of initial conditions, so that obseryables should not 
depend on the initial momentum spectrum proyided the 
energy density is fixed. Similarly, the momentum distri¬ 
bution in thermal equilibrium is uniyersal, therefore one 
expects obseryables to be independent of the differential 
cross section in this limit. 





FIG. 1. Left panel: VnifPr) as a function of l/ctot (thermal 
initial distribution, isotropic scattering cross section) in the 
rapidity interval \y\ < 0.5. The y-intercept is the ideal hydro- 
dynamic limit , while the slope corresponds to the viscous 
correction 5vn, as defined by Eq. ([T|). 



FIG. 2. Ideal hydrodynamic limit for the pr spectrum in the 
rapidity interval \y\ < 0.5. The spectra have been divided by 
dN/dy so that they are normalized to unity. 
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FIG. 3. Ideal hydrodynamic limit of anisotropic flow Top 
to bottom: n = 2, 3, 4. For all harmonics, we have chosen 
En ~ 0.2, SO that the linear response coefficients = Vnlsn 
are larger by a factor 5. 

Figure [2] displays the ideal hydrodynamic limit for the 
transverse momentum distribution in the rapidity inter¬ 
val \y\ < 0.5. As expected, it is independent of initial 
conditions and of the differential cross section. The dis¬ 
tribution is essentially exponential in pT, with a slight 
upward curvature which is a typical consequence of col¬ 
lective motion ini- The ideal hydrodynamic limit of 
anisotropic flow coefficients is displayed in Fig.|31 Vnipr) 
is linear at high px, as observed in previous numerical 
calculations fllll. This behavior is general in ideal hy¬ 
drodynamics [33| in the linear regime £„ 1. Inter¬ 

estingly, the magnitude of the linear response Knipr) = 
Vn{PT)l£n [ll] is approximately the same for n = 2,3, 4 . 
It would be interesting to reproduce these results using 
ideal hydrodynamics with the same ideal gas equation of 
state (constant sound velocity Cg = I/a/S) and a small 
freeze-out temperature. 


IV. VISCOUS CORRECTION 

We now present results for the first-order viscous cor¬ 
rection to observables, corresponding to the term Sf in 
Eq. m). We scale the viscous correction 6f by the ideal 
hydrodynamic limit Eq. ([1]) shows that Sf has 



FIG. 4. First-order relative viscous correction to the normal¬ 
ized transverse momentum distribution in the rapidity inter¬ 
val |i/| < 0.5. Dashed lines are power-law fits. 


the dimension of the cross section dtot ■ The relative vis¬ 
cous correction is 5///*^°^ divided by the total cross sec¬ 
tion tTtot- As a rule of thumb, viscous hydrodynamic 
applies if is significantly smaller than crtot in ab¬ 

solute value. 

Figure m displays Sf/f^^^ for the transverse momen¬ 
tum spectrum. Viscous effects result in a particle excess 
at large px, corresponding to an increase of the aver¬ 
age px, that is, a higher effective temperature. The rea¬ 
son is that viscosity decreases the longitudinal pressure, 
thereby reducing longitudinal cooling |4§|. 

Unlike the ideal hydrodynamic limit, the first-order 
viscous correction depends on the differential cross sec¬ 
tion. It is larger for smaller values of ttid- This is due to 
the fact that the scattering is forward peaked and less ef¬ 
ficient in thermalizing the system. For an isotropic cross 
section, the relative viscous correction is almost linear in 
Px' Sf oc px^'^^. This result is consistent with the 
results obtained in (4^ in the Chapman-Enskog approx¬ 
imation. 

Surprisingly, the viscous correction also depends on the 
initial momentum distribution. Two initial conditions 
with exactly the same energy-momentum tensor lead 
to different first-order viscous corrections to observables, 
at variance with usual viscous hydrodynamics [T^ . The 
reason is that the constant momentum distribution m 
is strongly out of equilibrium, therefore hydrodynamics 
does not apply at early times. Figured shows that with 
this constant distribution, Sf is smaller at high px than 
with a thermal initial distribution. This depletion at high 
Px can be understood as a memory of the initial condi¬ 
tions, where all particles initially have momenta below a 
threshold po- 

Finally, we study the viscous correction to anisotropic 
flow Vnipx)- Figure [5] displays —Svn/vn'^ as a function 
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FIG. 5. —Svn{pT) / Vn'’ as a function of transverse momentum 
Pt in the rapidity interval |i/| < 0.5. Open symbols refer to the 
case of isotropic cross section and full symbols to a forward- 
peaked cross section. The circles, triangles and squares refer 
to n = 2, 3 and 4 respectively. 

of transverse momentum pT- Viscous effects decrease 
anisotropic flow 0 , therefore the correction is shown 
with a minus sign. The viscous correction to Vn in¬ 
creases as a function of harmonic order n, as already 
observed in viscous hydrodynamic calculations [HI- At 
large transverse momentum, it scales approximately like 
the order n [s^, —Svn{pT)/v^'^ oc n. This dependence 
is weaker than the dependence reported in previous 
studies [dll-IH^. 

As expected, the viscous correction depends on the dif¬ 
ferential cross section. As for the spectrum, it is larger 
with a forward-peaked cross section {rriD = 0.3 GeV) 
than with an isotropic cross section. Results with the 
thermal distribution ([S]) and with the constant distribu¬ 
tion ([HI) (not shown) are consistent within error bars, 
thereby suggesting that the viscous correction to Vn is 
roughly independent of the initial momentum distribu¬ 
tion. 

With an isotropic cross section, the relative viscous 
correction saturates at high pr, while the usual quadratic 
freeze-out ansatz [l3l| predicts a linear increase. A mild 
increase with px is observed with a forward-peaked cross 
section, but it is still much slower than linear. Note 
that our viscous correction must be compared with the 
full viscous correction in hydrodynamics, which results 
in part from the hydrodynamic evolution, and in part 
from the freeze-out. With a quadratic freeze-out ansatz, 
however, the behavior of Vn at large px in hydrodynam¬ 
ics is dominated by the correction at freeze-out. Thus 
the quadratic freeze-out ansatz currently used in most 


viscous hydrodynamic calculations is not supported by 
microscopic transport calculations. 

V. CONCLUSIONS 

We have calculated the ideal hydrodynamic limit and 
the first-order viscous correction for the transverse mo¬ 
mentum {px) distribution transverse-momentum spec¬ 
tra and anisotropic flow Vn{px) by solving numerically a 
Boltzmann equation and studying the limit of large scat¬ 
tering cross section atot oo. The ideal hydrodynamic 
limit is found to be independent of microscopic details, as 
expected from the universality of thermodynamic behav¬ 
ior. The linear response coefficients Vnipx)/£n depends 
little on harmonic n in the ideal hydrodynamic limit. 

The first order viscous corrections to observables, on 
the other hand, are not universal. As expected [l^, they 
depend on the differential cross section. For all the dif¬ 
ferential cross sections investigated in this paper, we find 
that the relative viscous correction to anisotropic flow, 
Vn, does not increase significantly with px at large px- 
Our results suggest that first-order viscous corrections 
do not explain the decrease of at high px, at variance 
with common lore [H, and that a different mechanism, 
such as jet quenching [s^, is needed at high px- 

The viscous correction to increases linearly with n. 
The stron ger (rt^) dependence typically found in hydro¬ 
dynamics leads to negative values of V 4 and at 
large px [l4| . even for small viscosities. A weaker de¬ 
pendence on harmonic n is therefore likely to improve 
agreement with experimental data. 

Finally, our results clearly show that usual, second- 
order relativistic hydrodynamics does not fully capture 
the first-order viscous correction to observables. Two ini¬ 
tial conditions with exactly the same energy-momentum 
tensor, which would therefore yield the exact same hydro- 
dynamical flow, are found to yield different momentum 
spectra at the end of the evolution. Specifically, the first- 
order viscous correction is found to retain the memory of 
the initial condition. Within a strong coupling calcula¬ 
tion, a proper treatment of the underlying microscopic 
degrees of freedom leads to an evolution equation which 
is second-order rather than first-order [ 5 ^, so that the 
solution is not solely determined by the initial value of 
Our calculation provides an explicit illustration, 
within a weak-coupling calculation, that the initial value 
of does not solely determine the evolution. 
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